Generation of Squeezed States of Nanomechanical Resonators by Reservoir 

Engineering 
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An experimental demonstration of a non-classical state of a nanomechanical resonator is still an 
outstanding task. In this paper we show how the resonator can be cooled and driven into a squeezed 
state by a bichromatic microwave coupling to a charge qubit. The stationary oscillator state exhibits 
a reduced noise in one of the quadrature components by a factor of 0.5 - 0.2. These values are 
obtained for a 100 MHz resonator with a Q-value of 10^ to 10^ and for support temperatures of T si 
25 mK. We show that the coupling to the charge qubit can also be used to detect the squeezed state 
via measurements of the excited state population. Furthermore, by extending this measurement 
procedure a complete quantum state tomography of the resonator state can be performed. This 
provides a universal tool to detect a large variety of different states and to prove the quantum 
nature of a nanomechanical oscillator. 

PACS numbers: 85.85.+J, 85.35.Gv, 42.50.Dv 



I. INTRODUCTION 

With fabrication of nanomechanical resonators with 
fundamental frequencies from 100 MHz up to 1 GHz jlj, 
y, yf the demonstration of their quantum nature, in 
particular, the creation of non-classical states of these 
mesoscopic systems, have attracted a lot of interest !3, 
y, S S H- Apart from the fundamental interest in 
the study of these systems, nanomechanical resonators 
are also of great importance for technical applications. 
While micron-sized cantilevers are already used to per- 
form atomic force measurements, their noise properties 
which set the limits for the sensitivity of force detec- 
tion |3, [i3 can be improved by scaling them down to 
the nanometer regime. Residual thermal fluctuations can 
then be reduced to the quantum limit by active cooling 
schemes as proposed in [ll|, [l^, [l^l • 

A further reduction of the quantum fluctuations below 
the standard quantum limit can be achieved by squeezing 
the resonator mode. The idea to use squeezed states for 
measurements beyond the standard quantum limit pl4| 
appeared first in the context of the detection of gravita- 
tion waves. In principle it applies to any system where 
a weak classical force which has to be measured acts on 
a quantum-mechanical oscillator. The measurement con- 
sists of three stages. In the first stage one prepares the 
oscillator in a squeezed state, so that the dispersion of one 
of its quadratures is reduced below the quantum limit. 
Next, one allows the measured force to act on the oscilla- 
tor. At last one measures the squeezed quadrature, with- 
out touching the other one. To apply the squeezing ideas 
to the nanomechanical systems, all three stages have to 
be implemented. In particular, for the third stage, one 
has to design a setup in which only one of the quadratures 
is being measured. In this paper we restrict ourselves 
to the first stage, i.e., the preparation of the squeezed 
state of a nanomechanical oscillator. We show that by 
coupling the oscillator to a Cooper Pair Box (Josephson 



qubit) and irradiating the system by bichromatic, phase 
coherent microwaves one can "cool" the oscillator down 
to the squeezed state. 

We start by reminding the reader about the basics of 
squeezing and discuss two different ways to achieve it. 
Then we analyze the coupled CPB - oscillator system and 
show that the cooling into the squeezed state is feasible. 
At last we discuss possible ways to detect the squeezed 
state. 



II. SQUEEZED STATES AND RESERVOIR 
ENGINEERING 



For a harmonic oscillator with a Hamiltonian H = 
hva'a, where a and a' are the usual creation and anni- 
hilation operators, a general class of Gaussian minimum- 
uncertainty squeezed states is defined by [l5( 



\a,e)^D{a)S{e)\0) 



(1) 



Here D{a) = exp(Q!a^ — a* a) is a displacement opera- 



tor, and 5(e) =exp(- 



-.n 



) denotes the squeezing 



operator. We refer to the state |e) = \a = 0,e) as the 
squeezed vacuum state. The absolute value of the com- 
plex number e = re*^ is called the squeezing parameter. 
For these squeezed states the quadrature components, 
Ai^2, defined by a = (Ai + iA2)e'^/2 f^jf^^ ^j^g ^jj^gj., 

tainty relation AA1AA2 > |, where the variance of one 
component, AAi = e~^ /2 is reduced below the standard 
quantum limit of 1/2, whereas the noise in the other com- 
ponent is enhanced, AA2 = e'^^ /2. This property can be 
exploited to improve the sensitivity of measurements. It 
is important to note that this asymmetric distribution of 
the noise is stationary only in a frame rotating with the 
frequency i^ of the harmonic oscillator, since an initial 
squeezed state |e) evolves in time as lee"'^"^*). Therefore, 
the error ellipse rotates in phase space with the oscillator 




\9) 




FIG. 1: Representation of the squeezed state | e = re'*) 
in phase space spanned by the dimensionless position and 
momentum coordinates. The error ellipse indicates the re- 
duced/enhanced fluctuations for the quadrature components 

X1/X2. 



frequency, v, as shown in Fig.^ 

There are several ways to generate a squeezed state 
of a harmonic oscillator. A familiar method in quantum 
optics [l6j is to use a parametrically driven non-linear 
potential corresponding to the Hamiltonian 



H = hva)a - ihXia^^e-"^"* 



a^e*'^"*) 



(2) 



By going to a rotating frame, i.e. transforming the time- 
dependence away and assuming a parametric pump field 
frequency ujp — 2v, the Hamiltonian is simply given by 
Hi ~ ihX{a^ — a^^). Starting from the ground state of 
the harmonic oscillator, |0)i the time evolution with Hi 
produces the squeezed state \e — 2Xt)i. The application of 
this method for mechanical resonators has been proposed 
in Ref. Q, but the requirements, a sufficiently strong 
nonlinearity to overcome the losses, and an initial state 
close to the ground state, are not easily met. 

A second method, which we will elaborate on below, is 
to "engineer" an appropriate coupling to the environment 
such that a dissipative dynamics drives the harmonic os- 
cillator into a squeezed state, i.e. we "cool" the oscillator 
mode to a squeezed state. This reservoir engineering has 
been first proposed in Ref. [13, [l3 in the context of ion 
traps, and has been experimentally implemented in part 
by the Ion Trap Group at NIST in Boulder [Hi. This 
reservoir engineering can be achieved, for example, by 
coupling the oscillator to a dissipative two level system 
(TLS), where the form of the coupling determines the 
stationary state. 

The simplest (although trivial) example is provided by 
a reservoir which cools the oscillator to the ground state 
|0). We assume that the oscillator is coupled to a two- 
level system with ground and excited state \g), |e) ac- 
cording to the Hamiltonian Hi = g{aa^ + a)a^)^ where 
we use Pauli spin notation a+ = \e){g\ etc.. We further- 
more assume that the two-level system decays from the 
excited state to the ground state with a rate F. The time 
evolution in the interaction picture is then described by 



FIG. 2: Coherent and incoherent processes as described by 
master equation Q. The excitations of the qubit coher- 
ently decrease ((7i,"red sideband") and increase ((72, "blue side- 
band") the phonon number of the oscillator. The energy of 
the excite state is then dissipated with the decay rate of the 
TLS, F. These two processes relax the system into the de- 
sired state, p = \e,g){e,g\ while the coupling to the thermal 
phonon bath, 7 thermalizes the oscillator state. 



the master equation 



dp 
'dt 



-i[Hi,p] 



(2(T_/9cr+ — a+a^p — p(T+cr_) . (3) 



For long times the system thus evolves to the steady state 
p = |0)(0| (g) |g)(g| since \Q)\g) is a "dark state" of the 
system Hamiltonian, i.e. Hi\0)\g) — 0. 

For a general Hamiltonian of the form Hi — Fa+ + 
F^a^ where F is a function of a and a^ only, the dy- 
namics of master equation © "cools" the system into 
the state p = |V')(V'I ® |5)(5'I- The stationary oscillator 
state IV') is determined by F\^) = 0. A squeezed vacuum 
state, |e = re*^) obeys the relation 



(acosh(r) -t- a' sinh(r)e 



0, 



(4) 



and thus we choose F to be of the form [I3 

F = gia + g-ie-'^a^ , (5) 

where gi and (72 are related by r = atanh((72/5i)- For 
a single trapped ion driven by laser light and decaying 
via spontaneous emission such a Hamiltonian, Hi with F 
given in Eq. ((SJ can be constructed by applying two laser 
beams, one detuned to the "red sideband" (w = u>ge — v) 
and a weaker one detuned to the "blue sideband" (cj — 
^ge + ^) of the two-level transition frequency LOg^. This 
is illustrated in in Fig. |21 For a detailed explanation the 
reader is referred to the review article by Leibfried et 
al. p0|. 

In the following section we show how such a Hamil- 
tonian can be realized for a nanomechanical resonator 
coupled to a Cooper pair box, which plays the role of the 
dissipative two- level system [l^ . 

For a mechanical resonator the coupling of the oscilla- 
tor mode to the finite temperature phonon bath of the 



support has to be taken into account, leading to addi- 
tional contributions in the master equation. As we dis- 
cuss below, in the limit of high resonator frequencies 
where the rotating wave approximation (RWA) is valid 
the master equation is 



dp T 

— = -t[Hi,p\ + - [2a-pa+ - a+a-p- pa+a-) 

+ {Nb + 1)^ {2apa^ - a) ap - pa) a) 
H- -^B— (2a} pa — aa^ p — paa^^ . 



(6) 



The new terms in the master equation describe the heat- 
ing and dissipation of the resonator with rate 7, which 
relaxes to a thermal state with mean phonon occupation 



it. 



Nb- This will degrade the squeezing. 



III. THE MODEL 

We consider a nanomechanical resonator which is 
placed close to a Cooper pair box (CPB) as shown in 
Fig. 121 Similar systems have been considered in the 
context of cooling il3| and the generation of entangled 
states • The energy spectrum of the CPB alone is con- 
trolled by the gate voltage, Vg and the Josephson energy, 
Ej [2j. To obtain a coupling between the two systems 
the voltage T4 is applied on the resonator, leading to a po- 
sition dependent interaction via the capacitance Cx{x). 
The whole system is described by the Hamiltonian 



H = 



(Q-Q. 



2Cj 



\2 

Ej cos(0) + hva^a . 



(7) 



Here Qg — CgVg + CxVx is the total gate charge, Ce is 
the total capacitance of the island and v is the frequency 
of the fundamental flexural mode of the resonator. We 
decompose the voltages into a sum of a constant and 




FIG. 3: Setup: The mechanical resonator (thin, vertical bar) 
is biased by the voltage 14 to achieve a coupling to the Cooper 
pair box (CPB) via the capacitance Cx- The state of the CPB 
can be controlled independently by the gate voltage, Vg and 
the Josephson energy, Ej. 



a time dependent part, Vi = V^ + Vi{t) and therefore 
Qg = 2e(ng -|- ng{t)). For the desired interaction and to 
minimize relaxation (see below) the system is operated 
close to the degeneracy point, n^ = n -I- 1/2, ng{t) « 1 
and we reduce the CPB to an effective two level system 
with the basis states 



\g) = {\n) + \n+l))lV2, 
|e) = (|n)-|n+l))/x/2. 

The Hamiltonian simplifies to 

Ej 

H = + -TT'^z + 4:Ec{x)ng{t,x)(Tx 

+ Ec{x){l + 4ng{t, x)^) + hva)c 



(8) 



(9) 



Note that Ec{x) = e^/2CE as well as ng(t,x) depend on 
the resonator coordinate x via the capacitance Cx (x) . If 
we expand Cx up to the first order in x/d, where d is the 
distance between the resonator and the CPB we obtain 

E^ix) =. E, + E,^^ + Oi^), 

ng{t,x) = l.^(CgVgit)+Cyx{m-^))+OC^). 

Re-substituting these expansions into equation ((SJ and 
absorbing small shifts of the equilibrium position of the 
resonator in a redefinition of a and a) we obtain the fol- 
lowing contributions. 

H=+^o,^ + ^iCgVgit) + Cyxit))ax 



4^c /C° 
2e \Cs 



^CgVgit) - C^xim 



Cy^ 



d^ 



(10) 



The second term in Eq. H10() leads to direct excitations 
of the charge qubit. To avoid these excitations we choose 
the voltage signals such that they do not alter the gate 
charge, Qg{t)/2e = CgVg{t) + ClVx{t) w 0. For the sys- 
tems under consideration the expansion parameter, {x)/d 
is in the order of 10^^ and we can neglect all higher order 
terms in H. 

With these approximations the system Hamilto- 
nian (|10() reduces to 

H ^+^az- \{t){a + a^)ax + fiva''a. (11) 

Up to now the result is valid for arbitrary driving sig- 
nals. For the generation of squeezed states we choose the 
driving voltage to be of the form 



'{Ej-v)t\ f{Ej + v)t 

Vx{t) = Vi cos ( ^ ' ^ M + V2 cos [ \ 



(12) 



It consists of a part tuned to the red sideband and a 
part tuned to the blue sideband of the qubit transition 
frequency related by a fixed phase difference, 9. Finally 
we perform a transformation into the interaction picture 
with respect to Hq — Ej/2az + hva^a. Under the as- 
sumption |A(t)| ^ hv the RWA can be applied and we 
end up with 

H = h {gia + gze-'^a^) (7+ + h {g2e'^a + .gja^) a^ , (13) 



with the parameters 



-2E, 



' d 2e 



(14) 



where Xq = ^Jhj^mv is the extension of the resonator 
ground state. 

Discussion. For typical parameters values C^ « 2 x 
10"" F, Ec « 40 GHz, xo/d « 10"^ we obtain a value 
for the coupling strength of about g^ « 5 MHz for driving 
voltages still below 1 V. This value is also consistent with 
the approximations we have made (ng{t) <^ 1, RWA, . . . ), 
assuming a resonator with a fundamental frequency i/ > 
100 MHz. 

In practice, a perfect realization of the balance condi- 
tion Qg(t) = is impossible which leads to direct exci- 
tations of the charge qubit. However, an accuracy in the 
control of the voltages in the order of xo/d is sufficient to 
neglect this term since the applied voltages are detuned 
from the qubit transition frequency by v. 

Insufficient precision in the knowledge of the oscillator 
frequency as well as the qubit transition frequency leads 
to unavoidable detunings for the applied driving fields. 
Their effect is taken into account by adding the terms 
Sxa'a + Scq\e){e\ to Hamiltonian (|13|) . A detuning from 
the exact resonator frequency, Sx destroys perfect squeez- 
ing because the ideal state |e) is not an eigenstate of a^a. 
Since a measurement of the resonator frequency with a 
resolution of a few ppm can be achieved |3|, S^ is less 
than 1 kHz. The residual imperfection is a small effect 
compared to the influence of the finite Q- value and can 
therefore be neglected. The detuning from the charge 
qubit transition frequency, Scq is less crucial since it does 
not affect the steady state. For Scq < gi it only slightly 
changes the excitation probabilities of the qubit. A mea- 
surement with the required precision has been reported 
by Vion et al |2^. 

Damping. Apart from the unitary evolution given by 
H, the coupling to the environment provides the dissipa- 
tive part of the system dynamics. While a finite decay 
rate of the charge qubit is crucial for reservoir engineer- 
ing the damping of the resonator mode sets the limits 
of this method. Here we give a brief discussion of the 
dominant effects in our system due to the influence of 
the environment. 

The mechanisms of dissipation in superconducting 
qubits have not yet been fully investigated. The early ex- 
periments |23j| reported decoherence times of order sev- 
eral nano-seconds. This has been attributed to the ef- 
fect of the low frequency (1//) noise. In Ref. f22| it 



was demonstrated that this effect can be substantially 
reduced by operating at special symmetry points (de- 
generacy points). The decoherence time of 500 ns was 
achieved. Moreover, it became clear that a substantial 
part of the decoherence at such points is due to the energy 
relaxation (Ti in NMR) processes. Here, for simplicity, 
we assume that only the energy relaxation is important 
at the symmetry point. It is provided by the high fre- 
quency modes of the environment. One of the possible 
relaxation channels is the electro-magnetic environment 
in the external circuits which created fluctuations of the 
gate voltages. Taking into account these fluctuations by 
substituting Vi — > Vi + SVt we obtain the additional terms 
^§f-CiSViax due to this voltage noise. Assuming equal 
noise characteristics for SVg and SVx this implies a decay 
rate 



ca 



ft2 



CI 



(J2 

-^Svi+Ej/h). 



(15) 



For an external impedance Z{w) the noise spectrum 
of the voltage is given by S'y(aj) = 2ReZ{w)fihj[l — 
g-fic^/fesTj-i^ Because the temperatures reached with di- 
lution refrigerators are in the order of 10 — 50 mK which 
is much smaller than Ej /ks excitations of the qubit can 
be neglected and the decay rate simplifles to 



^x + C'g R Ej 



Ro h 



CI 



(16) 



where Rq = h/Ae^ is the resistance quantum. The decay 
rate can be adjusted by the gate capacitance, Cg and has 
typically values of 1-10 MHz 22). 

The dominant mechanism for the damping of the res- 
onator mode is the coupling to the phonon modes of the 
support. This leads to a finite decay rate 7 = i^ /Q^ where 
Q is the quality factor of the resonator. In contrast to 
the charge qubit the temperature of the environment is 
higher or comparable to the oscillator frequency. There- 
fore, the phonon modes at the resonator frequency have 
a non-zero occupation, Nb — [e'''^/'^^^ ~ 1]~^ and cause 
downward and upward transitions. For temperatures of 
10 — 50 mK the oscillator {v — 100 MHz) has an equilib- 
rium occupation number Nb ~ 2 — 10. 

Together with the Hamiltonian H13|l the decay rates F, 
7 and the bath occupation number Nb lead to a dissi- 
pative dynamics of the system described by the master 
equation 10. 



IV. RESULTS 

We are interested in the properties of the steady state 
solution of master equation ©. For a characterization 
of the stationary state we concentrate on the variance of 
the Xi quadrature component, (AA'i)^^ = {Xf) — (A'l)^ 
where the average is taken with respect to the steady 
state density matrix, ps. We compare the variance AATi 



with the zero point fluctuations, (AXi)q 
define the ratio 



7^ = 



(AX 



ijo 



]X^\0) and 



(17) 



as a measure for the degree of squeezing. 

To study the effect of the different terms in Eq. © it is 
convenient to look at the master equation in the squeezed 
frame, i.e. we perform the unitary transformation, U — 
S{e) where the value of e = re*^ is chosen according to 
Eq. lO. For the transformed density operator, p = W pU 
we obtain the equation 



dp .,., 



a^a' 



P\ 



+ 2 (2(^-P^+ - ^+<^-P - P<^+<^-) 



+ ^{N+l){'2apa) 



2 

7 



a* ap - 



pa' a) 



+ —N [2a' pa — aa' p — paa*^ 



(18) 



M {2apa — aap — paa) 



- ]-M* (2a^pa^ - 
where g — |.gi|/cosh(r) and 



a' a' p 



pa* a'^ 



TV = (/Vb + 1) sinh2(r) + N^ cosh2(r) 
M = (2A^B + l)e'^ cosh(r) sinh(r) . 

The master equation can be written as the sum of three 
Liouville operators 



jiP ^ (^^ 



Cr 



C^)p. 



(19) 



In the squeezed frame the first two terms correspond to 
the master equation that is known from sideband cool- 
ing in ion traps (Eq. (O). Excitations of the qubit on 
the red sideband followed by a spontaneous decay succes- 
sively reduce the phonon number and drive the oscillator 
into the ground state which corresponds to a squeezed 
vacuum state in the original frame, |e) =C/|0). The third 
contribution, £-y describes the coupling of the oscillator 
mode to a squeezed reservoir [23 . The degree of squeez- 
ing of the reservoir is maximal since the parameters N 
and M are related hy \M\'^ ^ N {N + \) . The different 
processes of the system dynamics in the squeezed frame 
are summarized in Fig. ^ Note that N and M grow ex- 
ponentially with the squeezing parameter while the Rabi 
frequency, g decreases exponentially. Therefore, even for 
a weak coupling, 7 the effects of the environment become 
essential as soon as r ^ 1. 

Weak coupling. We first consider the limit where the 
decay of the charge qubit is much faster than the rest of 
the system dynamics, F ^ g{r),jN{r). In this regime 
the excited state can be adiabatically eliminated by treat- 
ing the coupling term, Cg in second order perturbation 



l.9> 




-^^i^o/(7V+l) 



n = 



n=l n=2 



FIG. 4: Coherent and incoherent processes in the squeezed 
frame as described by master equation 1181 . In this picture 
the qubit is only excited on the red sideband {g). Followed by 
a spontaneous decay (F) these transitions cool the oscillator to 
its ground state. The cooling is compensated by the squeezed 
heat bath which apart from the enhanced heating rate (7A'') 
also induces coherences between the states |7i) and |n ± 2) 
(7M). 



theory. After tracing over the qubit degrees we obtain a 
new master equation for the oscillator density operator, 
Px ■ The excitations of the charge qubit on the red side- 
band provide an additional cooling rate of 4g^/F for the 
resonator. The master equation for the oscillator alone 
can be solved by transforming it into a partial differen- 
tial equation for the Wigner function. The details of this 
calculations are given in Appendix ^ In this limit we 
obtain the result 



n = 



2s!. 

r 



-^^ + ^{2Nb + 1) 



\ Mi I 2 

\ r ^ 2 



(20) 



For small r and 7 this is close to the ideal behavior TZ = 
e"*" wh ile for high values of r it saturates at the value of 
TZ = yJ{2NB + 1) which corresponds to the unperturbed 
thermal state. 

Strong coupling. In the strongly driven regime, g S> 
F, 77V, the system performs rapid oscillations between 
the states \n,g) and \n — l,e) on a timescale which is 
much faster than the incoherent processes. It is there- 
fore convenient to work in the basis of dressed states 
|n,±) = ^(|n,0)±|n-l,l)),|0) = |0,0). Since these are 



V2^ 



the eigenstatcs of the Hamiltonian H = g{a(j^ + a* a^) 
the steady state density operator of Eq. IjlSI) becomes 
diagonal in that basis as ^ ^- 00. By neglecting the off- 
diagonal terms the density operator can be approximated 



p = Po|6)(0|+^p+|n,+)(n,+|+p-|n,-)(n,-|. (21) 

Because the Liouville operators Ct and C^ do not dis- 
criminate between the states \n, ±) we define the joint 



probabilities Pn 



- p„ . With this ansatz the master 



(22) 



equation H18|l reduces to the rate equation 

■Pn= + T+in - l)pn-l + r_(n + l)pn+l 

-[r+(n)+r_(n)]p„, 
with the heating and coohng rates 

T_H = (r + 7(iV + l)(2n - l))/2, r_(o) - , 

T+{n)=-fN{2n+l)/2, T+{Q) = -fN . 

(23) 

In the stationary state the occupation numbers are deter- 
mined by the detailed balance condition T_(n+l)p„+i — 
T^(n)pn and an analytic expression for the mean occu- 
pation number is given by 



{a^a)p = 



jN (2 2Fi[|, 2, 1 -(- a, z] - 2^1 [1, f, f + «, z]) 



r + -f{N+l) + jN2Fi[l 



3 3 

2' 2 



■a,z] 



. (24) 

2i^i denotes the hypergeometric function depending on 
the parameter a — T /{2'^{N + 1)) and is evaluated at 
the argument z = N/{N + 1). Since {a?)~p, {a'^^)p -> 
in the limit of strong coupling, a transformation back to 
the original frame simply gives 



7^ = e" 



2(ata)^ 



1. 



(25) 



Perturbation theory. The interesting regime where the 
final state of the resonator is squeezed, TZ < 1, obviously 
requires 7iV ^ g, T. With this restriction a solution for 
arbitrary parameters g and T can be found by taking the 
ideal solution, p = |0, g) (0, g\ and treat the corrections of 
Cy in first order perturbation theory. The details of the 
calculations are listed in Appendix ^ and the result of 
this approach is 



7^ 



^-*(f 



r 



(26) 



While this expression gives the correct interpolation be- 
tween the weak and the strong coupling limit it is only 
valid for (a^'a)p <C 1. A rough estimation shows that this 
is still true up to the minimum, 7?,,nin = min{7?,(r), r > 0} 
in the case of 51 < F while expression (|26|l gives rather 
poor results for 7?.min in the case of gi > F. 

Numerical results. For numerical calculations the mas- 
ter equation H18I) is evaluated in the number basis. In the 
squeezed frame the solution is close to the ground state so 
a relatively small number of matrix elements is sufficient 
to describe the exact state. Fig.|21shows numerically cal- 
culated values of 7?.,nin for various parameter values for 
gi, F and 7. The results show that the noise, AXi can be 
reduced to half of the standard quantum limit for damp- 
ing rates 7 > 0.02. This corresponds to a Q-factor of 
5000 in the case of a 100 MHz resonator. For Q = 10^ 
a reduction by a factor of 5 is possible, still assuming a 
"hot" environment of about T w 20 — 30 mK. Obviously, 
higher oscillator frequencies or lower temperatures would 
improve the results even further. 
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FIG. 5: Minimum values for the noise reduction, 7?.min for var- 
ious parameter values, gi, F, 7 (all in MHz) and Nb = 5. Fig- 
ure a) shows the dependence of 7?.min on the driving strength, 
g\ for a fixed F = 2, while b) shows the dependence on the 
decay rate of the charge qubit, F for fixed (?2 = 5. The difi'er- 
ent values for 7 in a) and b) are 0.001 (solid), 0.005 (dashed), 
0.01 (dotted) 0.02 (dashed- dotted). In Figure c) 7?.min is plot- 
ted as a function of the phonon damping rate, 7 for g\ = ^ 
and F = 2 (solid), F = 5 (dashed), F = 8 (dotted). 



V. DETECTION OF A SQUEEZED STATE 

Below we discuss schemes for detecting non-classical 
states, in particular squeezed states, of a nanomechan- 
ical system. The displacement measurements based on 
laser interferometry, which are used for mechanical res- 
onators with a length in the order of a hundred mi- 
crons, cannot be applied in the nanometer regime. Al- 
ternatively displacement detectors based on a single elec- 
tron transistor (SET) [23 have been considered 2^ 2iM 
and were recently used by the groups of A. Cleland [23 
and K. Schwab |29| to measure the fluctuation spectrum 
of nanomechanical resonators. While in current experi- 
ments the displacement sensitivity is still limited by the 
amplifier noise, the quantum limit is determined by the 
back action of the current shot noise on the resonator. 
By increasing the signal amplification of the detector, 
which is necessary to observe the reduced fluctuations of 



a squeezed state also this back action is enhanced. A 
quantum mechanical analysis of the properties of SET- 
based displacement detectors is presented in J2^|23- Us- 
ing the results of the analysis done by Mozyrski et al. |2y| 
the charge fluctuations on the SET island would destroy 
the squeezed state, especially in the experimentally at- 
tractive sequential tunneling regime. 

In quantum optics with trapped ions and in Cavity 
QED experiments, information about the oscillator state 
is often obtained via a coupling to a two level atom. Ef- 
ficient readout techniques for the state of an atom can 
be used to measure properties of the not easily accessible 
oscillator mode. Since CPB and other TLS are currently 
developed for quantum computation |2j |22, ISJ , which 
in particular implies a read out of the qubit represented 
by e.g. the charge states of the CPB, this "measurement 
toolbox" is being developed in mesoscopic physics. Mo- 
tivated by this, we discuss the detection of the resonator 
state via the readout of the charge qubit. We concen- 
trate on two detection methods that can be performed 
with the setup shown in Fig. |2| extended by a measuring 
device for the state of the CPB. 



A. Dark Resonance 

A simple way to verify the generation of a squeezed 
oscillator state is to look at the excited state population, 
Pe of the charge qubit. Since the squeezed state is a dark 
state of the system Hamiltonian Ijl^f) the qubit excita- 
tions are significantly suppressed even in the presence of 
the driving fields. By varying the detuning 5x a dark res- 
onance becomes visible bX 5x = Q [iTij which corresponds 
to the generation of a squeezed state. Fig.Elshows the ex- 
pected correlations between the degree of squeezing and 
the steady state excitations of the qubit as a function of 
the detuning 5x- 

In the presence of a finite 7 the population pe retains 
at a value of about jN/T ^ 1 (see Appendix IBl) . In the 
regime of strong coupling, (71,52 3> F this is clearly dis- 
tinguishable from the value Pe ~ 1/2 as expected for, e.g. 
a thermal state. For weak driving fields and a thermal 
oscillator state we expect a excited state population of 
Pe w {n)gl/r^ + ((^) + Ijffl/r^- Therefore the condition 
5I/F > 7A^ to distinguish the squeezed state from a low 
temperature thermal state. 



B. Occupation Numbers 

According to its definition (Eq. (^) a characteristic 
property of the (ideal) squeezed state, |e) is that only the 
even number states are populated. The measurement of 
the resonator populations via the Stark shift of the qubit 
resonance frequency was proposed in Ref. [3 while in 
Ref. [331 it is suggested to utilize the anharmonic coupling 
between bending modes for a Fock state readout. Here 
we follow a different line |33| and use the linear coupling 




FIG. 6: Correspondence between the excited state popula- 
tion, Pe {solid line) and the degree of squeezing, 7?. {dashed 
line) as a function of the detuning, Sx- The values are plotted 
for the parameter values 7 = 0.01 MHz, r = 0.7 and Nb = 5 
for a) gi = 5 MHz, F = 2 MHz {strong coupling) and b) 
gi = 1.5 MHz, r = 5 MHz {weak coupling). 



to the TLS to determine the occupation numbers, Q„ — 
{n\px\n). The basic idea is as follows. Suppose we start 
from an initial density operator p(0) — Px ® \9){9\ ^^d 
switch on a Jaynes-Cummings coupling. Hi = hg{a<7^ + 
a^a-) between the charge qubit and the resonator. The 
evolution of the qubit polarization is then given by 



{a,){t) ^ -J2QnCOs{2nnt) . 



(27) 



Due to the different Rabi frequencies Qn = g\/n the val- 
ues of Qn can be extracted from the Fourier transform 
of this signal. 

In our system the required coupling of the oscillator 
to the qubit is realized by Hamiltonian (|13|l with gi = g 
and 92 = 0. In contrast to the ideal situation of Eq. H27|) 
the decay of the charge qubit and thermalization of the 
resonator mode lead to modifications of the signal and 
restrict the applicability of this method. 

Obviously, a necessary condition to resolve the oscil- 
lations of the qubit polarization is the strong coupling 
regime g ^ T,j{Nb + 1). An approximate time evo- 
lution of the system can be obtained by the following 
considerations. Starting from a pure state \n,g) the sys- 
tem will oscillate between this state and \n — l,e) with 
a frequency 2ri„, where J7„ = g\/n- During this oscil- 
lation it decays into neighboring number states with a 
rate i?„ = r+ {n) + T_ (n) as defined in Eq. ^^ . Since 
all other states, |m ^ n, g) and \m ^ n — l,e) are popu- 
lated gradually their oscillations wash out and with the 



exception of the ground state they give no contribution 
for az{t) ■ Therefore, for the initial state \n,g) we obtain 



{o-z){t) = -po{n,t) - cos(2fi„i)e" 



-R,,t 



(28) 



where poit) is the population which accumulates in the 
ground state which has the form 



PoKO = l-e-"/^Ep 



fc=0 



kl 



(29) 



in the limit F ^ 7(-/Vb + 1). The exact time dependence 
for 7 > is not important since the changes of po("'j t) are 
slow compared to the Rabi oscillations. For an arbitrary 
initial state with occupation numbers Qn the polarization 
of the qubit is given by 

{a,){t) = -Y,Qn{po{n,t)+cos{2nr.t)e-'^"*) . (30) 

n 

The extraction of the occupation probabilities from a 
given function {fjz){t) requires the resolution of the indi- 
vidual Lorentzian peaks in the Fourier transform of this 
signal. The n— th peak can be resolved if the condition 
Rn + Rn+i < gl \fn is true. Therefore, a lower bound 
for the maximum occupation probability which we can 
determine with this method is given by the solution of 
the equation 

(F + 2^NB)n\l:t^ + 27(27Vb + 1)^11 = g ■ (31) 

A first order approximation, valid for ^g^ /V"^ <C 1 gives 






A-f{2NB + l)g^ 
F3 + 67(2iVB + l).g2 



(32) 



For the parameter values g = 5 F, Nb = 5 and 7 = 
0.001 — 0.01 F we obtain Umax ~ 8— 17, which is sufficient 
to detect slightly squeezed or low number states. 

Fig. [3 shows the time evolution of {<Jz){t) and the 
extracted occupation probabilities, Qn for a non-ideal 
squeezed state. A clear distinction between the squeezed 
state and a thermal state with the same mean occupation 
number is possible. 

Although this method provides a complete determi- 
nation of the number state populations for states close 
to the ground state, it is not possible to distinguish be- 
tween coherent superpositions and mixed states. There- 
fore, in the next section, we discuss the implementation 
of quantum state tomography to obtain the full informa- 
tion about the resonator's density matrix. 



VI. QUANTUM STATE TOMOGRAPHY 

The ultimate determination of an arbitrary quantum 
state is the measurement of the complete density oper- 
ator. The procedure of estimating the density matrix 
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FIG. 7: Measurement of the oscillator occupation numbers, 
Qn- A squeezed state is generated as explained in the previous 
sections with (all in MHz) g\ = 5, 32 = 3.5 (r = 0.88), F = 2, 
7 = 0.01 and Nb = 5 MHz. For the detection we assume 
the same values but g = gi =8 and (72 = 0. Figure a) shows 
the exact time evolution of (crz)(t) (solid) which is compared 
to the expected behavior (dashed) as given in Eq. 1281 . The 
difference for larger times is due to the deviations of the exact 
time dependence of po(n,t) from Eq. 12911 for 7 > 0. In Fig- 
ure b) the extracted occupation probabilities (dark gray) are 
compared with the exact values (hght gray). The inset shows 
the estimated number distribution for a ideal squeezed state 
(top) and a thermal state with the same mean occupation 
number, (n) = 1 (bottom). 



by repeated measurements on the same initial state is 
called quantum tomography .34] . For a harmonic oscil- 
lator the Wigner function of a state, a quasi probability 
distribution in phase space J2J, l35| . contains the same 
information as the density matrix. The implementation 
of a method to reconstruct the Wigner function for a 
nanomechanical resonator provides an universal tool to 
detect and characterize various non-classical states, and 
therefore a tool to clearly demonstrate the quantum na- 
ture of this still macroscopic system. 

In quantum optics methods for state tomography are 
well known |35Ll36ll37l | and have been successfully imple- 
mented to detect non-classical states of a cavity field ,33| 
or the motion of a trapped ion [33 ■ We consider the 
method discussed in Ref. "331 and show that it is appro- 
priate for the implementation in nano scale mechanical 
systems. 

To determine the Wigner function at a certain point 
in phase space, a = x + ip we start from the identity |35j 






(33) 



n=0 



where Qn{ct) — {nlD^^ {a)pxD{a)\n) are the occupation 
probabiUties of the displaced density operator. The val- 
ues of Qn{ct) are measured as discussed in the previous 
section. 

A complete state tomography consists of the following 
steps. For each point in phase space (x,p) we displace 
the original density operator by a = —{x + ip). Then 
the occupation numbers of the oscillator are determined 
by measuring the time evolution of the qubit polariza- 
tion, ((Tz)(t). In the end we obtain the measured value 
of the Wigner function, VF(a) by summing up the Qn(a) 
according to Eq. (|^ . 

In the following we discuss the limits for the implemen- 
tation of the individual steps for a system consisting of a 
nanomechanical resonator coupled to a charge qubit. 

Applying the displacement operator. The first step of 
the procedure described above requires the shift of the os- 
cillator by the complex amplitude a. In the setup shown 
in Fig. 131 this displacement can be achieved either by ap- 
plying an additional voltage to a lead opposite the CPB 
or by exploiting the existing coupling to the charge qubit, 
X{t){a + a^)ax- 

In the first case a voltage drop Vd over a capacitance 
Cd{x) formed by the lead and the resonator generates 
the driving Hamiltoninan, Hd = X{t){a + a') with A — 
^CdV^Xo/d. The evolution of the oscillator under H — 



hva) c 



with 



Hd is m 



|^(i)) - b{a{t))e- 



'IVXO)), 



a{t) 



dt' e"'^'-'"^\{t') 



(34) 



(35) 



Because A ^^ 1 GHz ^ ;/, a short constant voltage pulse 
is sufficient to obtain displacements of |a| < 10. 

If the coupling to the CPB is used for the displacement, 
we first transfer the ground state of the qubit |0) into 
one of the eigenstates of ax by adiabatically changing the 
CPB parameters Vi* and Ej. Since the coupling strength, 
A « 5-10 MHz is much lower than the oscillator frequency 
a radio frequency pulse A(i) ^ cos[vt) has to be applied 
to achieve shifts in the order of \a\ > 1. 

Especially in the second case the Wigner function of 
the oscillator is modified during the displacement due to 
thermalization with the phonon bath. The diffusive dy- 
namics of a driven oscillator and therefore the resulting 
errors can be calculated exactly (see Appendix lU)). If we 
assume an initial Gaussian distribution and a displace- 
ment time At = |q;|/|A| ^ I/7 we obtain the relative 
errors for the widths 



eiKp) 



|A| 



1 



Nb 



2A2,^ 



(36) 



The damping also modifies the displacement amplitude. 
Together with the deviation caused by a small detuning 
in the driving frequency lu = 1/ + 5x, we obtain a relative 
error 



I IN |a| f\Sx\ ^ 7 



(37) 



For the parameter values considered in this paper 
e{\a\) ~10-3. 

Occupation numbers. The measurement of the proba- 
bilities, Q„(a) is done as discussed in the previous sec- 
tion. 

For an estimation of the error due to truncation we 
consider an initially density matrix of a pure number 
state, p = |?7i,)(?n| with m < rimax- The application of 
the displacement operator shifts some part of the wave- 
function out of the detectable subspace {|0), . . . , |nmax)}- 
This leads to an absolute error in the estimated Wigner 
function iy(a) of 

2 °° 
e{m,a)^\W{a)-W{a)\<- V \{k\D(-a)\m)\'^ . 

IT ^ — ^ 

fc — rimax 

(38) 
The matrix elements of the displacement operator are 
given by 



k\b{- 



fc! 



i"r 



[Ll 



where L^™ are the generalized Laguerre polynomials. 
The sum '^^{'m\p\m) £{m,a) <C 1 provides an upper 
bound for the total truncation error and rough estimation 
whether the method is applicable or not. In practice the 
determination of the Qn is done by fitting the actual 
signal by minimizing the deviations in the least square 
sense. For a linear fit the standard deviation of the Q„ 
can be written as 



c{Qn) = CnCFD ■ 



(39) 



The error in the measured data point, ao includes the 
error from the measurement of the qubit polarization as 
well as deviations of the system evolution from expected 
evolution as given in Eq. (|28|l . The coefficients, c„ de- 
pend on t he p arameters of the system and the fitting 
procedure po]. While the c„ increase on a scale set by 
'T-max it is still possible to determine the Q„ beyond rimax 
in the expense of accuracy. For the example given below 
the c„ are shown in Fig. |H1 

Example. To summarize the considerations made in 
this section we discuss the generation of a squeezed state 
and the reconstruction of its Wigner function for a spe- 
cific example in some detail. We consider a nanomechan- 
ical resonator with a fundamental frequency of 100 MHz 
and a Q-value of 2 x 10* cooled to a temperature of 25 
mK. For the coupling to the CPB we assume the val- 
ues 51 = 5 MHz, 32 = 3 MHz (r = 0.7) and a decay 
rate P = 2 Afflz. Using the results of the last section. 
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FIG. 8: Coefficients, Cn which determine the error of the esti- 
mated occupation numbers according to Eq. 1391 . For a large 
number of measured data points, Nd the Cn are proportional 
to l/\/ Nd- See text for details. 



this drives the resonator into a state with a noise reduc- 
tion TZ ~ 0.5 in one quadrature components. After the 
resonator has reached its steady state we switch off the 
coupling and wait for a short time t ^ l/F to reduce 
the excited state population to less than 1%. According 
to Eq. (|36|l this delay leads to a broadening in the Xi 
direction of about 17%, so we end up with TZ ~ 0.6. The 
exact Wigner function for this state is well located in a 
the phase space region |a| < 2. 

For the displacement of the oscillator we use the cou- 
pling to the CPB as described above with A = 8 MHz. 
Eq. ||2ni) predicts an error of ^(A^) « |a|x2% 

For the determination of the occupation numbers, 
(5ri(a) we apply a red sideband signal with g — S MHz. 
For this value we obtain n^ax ~ 8. The estimation 
'^^{m\p\m) e{m,a) < 0.15 for the phase space region 
|a| < 2 shows that an accurate reconstruction of the 
Wigner function is possible. To resolve oscillations of 
frequency g^rimax over a time much longer than the char- 
acteristic decay time the number of measured data point, 
No has to fulfill No ^ 2g.y/nmax/(7''F). In our example 
we choose a measurement time T = 3 F^^ and Nd — 150. 
For these values the coefficients c„ in Eq. (|39|l are plotted 
in Fig.lHl We suppose that the qubit polarization can be 
measured with an accuracy of 0.02 and that the values 
for Q0...12 (we choose rimax = 12) are estimated. Adding 
the errors of the occupation numbers, the error due to 
truncation and the error from the displacement we ex- 
pect an accuracy for the reconstructed Wigner function 
oi\W{a)~Wia)\ <0.05. 

The results of a numerical simulation of the generation 
of the squeezed state and the reconstruction of its Wigner 
function is shown in Fig. |5| 

In Appendix ^ we also briefly discuss a different 
method for a state to mog raphy as proposed by Lutter- 
bach and Davidovich J33. While this method is exper- 
imentally more attractive since it requires less measure- 
ments, a stronger coupling and longer decoherence times 
of the qubit are required for its implementation. 




FIG. 9: Results of the numerical simulation of the state to- 
mography for a squeezed state as described in Section IVII 
Figure a) shows the reconstructed Wigner function, W{a). In 
Figure b) a cut of W{a) (solid line) in x direction (Im(Qf) — 0) 
is compared with the corresponding values of the real Wigner 
function, IV(a) (dashed line) and the Wigner function of the 
oscillator ground state (dotted line). 



VII. CONCLUSION 

In this paper we showed that by reservoir engineer- 
ing the fundamental mode of a nanomechanical resonator 
can be driven into a squeezed state. The stationary state 
exhibits noise reduction in one of the quadrature com- 
ponents by a factor of7?.«0.5-0.2 below the standard 
quantum limit. For a 100 MHz resonator these values 
are obtained for Q-values in the order of Q = 10^ — 10^ 
and standard dilution refrigerator temperatures of T « 
30 mK. The detection of the squeezed state can be done 
within the same setup as used for its generation by mea- 
suring the excitation probability of the charge qubit. Fur- 
thermore this measurement procedure can be extended 
to obtain a complete reconstruction of the oscillator's 
Wigner function. This tool provides an universal detec- 
tion method for non-classical behavior of the resonator. 
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undo this transformation and finally obtain the Wigner 
function in the original frame, 



APPENDIX A: ADIABATIC ELIMINATION 



1 



Wix,p) = AAexp ( --(x,p)S"^(x,p)^ ) , (A5) 



We consider master equation H19|) in the limit F ^ 
g, "fN. To zero order in g/F, 7/F the equation reduces to 
dtp = CtP and a steady state solution is located within 
the subspace S — {px ® \g){9\}- The term C^ leads to 
a slow dynamic within this subspace while Cg couples S 
to its complement. Since Crp — OiV), \fp ^ S we can 
project the master equation (|18|l on the subspace S and 
treat the coupling of Cg in perturbation theory. After 
tracing over the charge qubit states we obtain 



where the variance matrix E is 



dpx 
dt 



C^Px ~ Trcgj/Ig^r ^ '^g P^® l5> (sll ■ (Al) 



An evaluation of this expression leads to a master equa- 
tion for the reduced density operator of the resonator, px 
given by 

— ^ = ( — + 2 (^ + 1) ) i'^O'Pxo' - a'apx - Pxo'a) 

+ -N {2a}pxa - aa) px ~ Pxaa^) 

T 

- -M (2apxa - aapx ~ Pxaa) 

- -M* (2a^pxa^ - a^a^ Px - Pxa^a^) ■ 

(A2) 

Because the master equation now only contains the op- 
erators a and a^ it can be transformed into a partial 
differential equation for the Wigner function JI3, 



d^ 



^ + i(^^+l)yaaaa* 



W{a) (A3) 



.1 ( M— + M*—— 1 W(a) 



By introducing the new variables i, p defined by a = 
(i -I- ip)e'*/^, Eq. IJA3|I separates into two independent 
Fokker- Planck equations and their steady state solution 
is 



J.-1 ^ rpT 



(T. ;.,)-. (-) 



with 



r^_( e^ W cos(0/2) sin(0/2)\ 
' e-' {- sin(6'/2) cos(e'/2) J 



The variances of the two quadrature components, X1.2 
are given by the widths of the Gaussian function 

W{Xi,X2), 



(A^.)^4(^e^^ + i 



(2...1))/(|^.i) 



(A7) 



APPENDIX B: PERTURBATIVE SOLUTION 

For a weak coupling of the oscillator to the phonon 
bath we write the master equation l|18(l as 



dp 



= Cq{p) + C^{p) 



(Bl) 



where Co = Cg + Cr contains the terms proportional to g 
and F with the steady state solution po = \0) {0\ ^ \g) {g\. 
In the limit of weak dissipation, 7(iV + 1) <C 5,F we 
write the density operator as p — po + pi, and calculate 
the corrections, pi up to first order in 7. 



dpi 
dt 



-Coipi) + CyiJJo) 



F, 



■i[g{aa+ -|-aV_) - i-|e)(e|, pi ] + Fcr_picr+ 

-7iV(|l)(l|-|0)(0|)+7M|0)(2H-7Ar|2)(0|. 

(B2) 



W{x,p) = A/'exp 



2g 



""^ l\ f i"^ P^ 



F 



Dx Dp 



(A4) 



where A/" is a normalization constant and 

Because these calculations have been performed in the 
squeezed frame we need another change of variables to 



The steady state solution can be calculated e.g. by eval- 
uating this equation in the number state basis. Using 
the notation P„ = {g,n + l\pi\e,n) — {e,n\pi\g,n + 1), 
En — {e,n\pi\e,n) and G„ = (g, njpijg, n), we obtain the 
set of coupled equations 



-ig^/^l + TPn - F£'„ 



0, 



-T/2Pn-i2gV^^Tl{E.n-Gn+i) = 0, 

ig^Pn-l+TEn+-iNSn,l''lN5nfi = 0. 
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Because all matrix elements for n > 1 are zero we find 
a very simple solution for the mean occupation number 
and the excited state population 



{n)^ = G,+E,=jN[- + — 



1 r 

f + 4^ 



{\e){e\)^^Eo^^ 



(B3) 



Another set of equations for the matrix elements between 
the states {n\ and |n+2) gives two non-zero contributions 
for {g,0\pi\g,2) and {g,2\pi\g,0) which lead to 



{a^')p^^{M + M*). 



(B4) 



The variance of the Xi quadrature component in the orig- 
inal frame is just given by 



iAX,r^'-^{{a'), + {a^'] 



p + 2{h)p + l) 



(B5) 



This leads to the result of Eq. if^ . 



APPENDIX C: DISSIPATIVE, DRIVEN 
OSCILLATOR 



We consider a harmonic oscillator, H — hva^ a driven 
by the linear term Hci = Xe^^^a + A*e^''^*a'f which is 
weakly coupled to a reservoir. The master equation in 
the interaction picture is given by 

— — —i[Xa + Aa' , p\ H yZapa' —a' ap— pa' aj 

I T^B /„ t t n 

-\ lyZa'pa — aa'p—paa'j . 

(CI) 

This equation can be transformed into a Fokker-Planck 
equation for the Wigner function [iJl- For the coordi- 
nates X — {a + a*)/2 and p = (a — a*)/2i we obtain 



dW 
~dt' 



R«(A)| + MA)|:' 



W 



d 



d 



dx dp 



Nb 



dx"^ dp^ 



W. 

(C2) 



A general solution for an initial distribution, Wi{a) is 
given by 



W{a,t)-- 



1 



27rCT2(i) 



d^Qfo Wi[aQ) exp 



I — ii 12 

\a — OfQe 2 1^ 



2a2(i) 



(C3) 
with 2cr2(i) = [Nb + \){l - e"^*). For the special case 
where the the initial distribution is a Gaussian function 



centered at the origin the time dependent solution can 
be evaluated as 



W{x,p,t) =7V(i)exp 



{x-cAt)f {p-cp{t)f 



2Ax^{t) 2Ap2(i) 



(C4) 

with a normalization factor, JV{t) and the time depen- 
dent parameters 



Ax^ft) = Ax^(0)e-^*+ ^^^/^ 



(1-e-^*), 



A/(i) = A/(0)e" 



-7t 



2iVij + l 



4 -d-^"^')' (C5) 
c,(f) = 2Im(A)(l-e-^*/2)/^^ 

Cp(t)=2Re(A)(l-e-^*/2)/7. 



APPENDIX D: FAST TOMOGRAPHY 

The reconstruction of the Wigner function as described 
in Section IVII requires a lot of measurements since it re- 
lies on the whole time evolution of the qubit polarization. 
A direct measurement of the Winger function form a sin- 
gle value of (uz) has been proposed by Lutterbach and 
Davidovich .373 foi' Cavity QED and ion traps. Here we 
give a brief summary of this procedure which relies on 
the identity 



W{a) = -Tr{D^(a)pxD(a)e''^^*^} , 



(Dl) 



Again the measurements on the qubit are used to deduce 
the expectation value of the parity operator, e"^" ° . To 
do so we assume that a certain evolution U can be applied 
to the system with the following properties: 

C^(l5>l^» - \9)\n) if n is odd, 
U{\g)\n)) = \e)\n) if n is even . (D2) 

Starting from an initial density operator p(0) — Px ® 
|0)(0| the Wigner function at the point a can then be re- 
constructed in three steps: First, apply the displacement 
operator, D{—a) — D^{a). Second, let the system evolve 
according to Eq. (|D2|I and third, measure the polariza- 
tion of the charge qubit. Then 

{a,) = Tr,+cQ{UDt(a)p(0)D(a)UVj 
= Tr,+cQ{D^(a)p(0)D(a)UV,U} 



Tr,{Dt(a)p,D(a)e'^ 



^} = -W(«), 



Due to the special properties of the evolution operator 
[/, the summation of Eq. H33|l is performed by the sys- 
tem itself and therefore only one value of (cr^) already 
determines W{a). 

The crucial point for the implementation of this 
method in mesoscopic systems is the construction of the 
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time evolution U as given in Eq. (ID2|I . It can be achieved 
with a Hamiltonian of the form ,37j 

H = -{Ej + hA{t)a^a)<T^+hiya'^a. (D3) 

Then a Ramsey interferonietry produces the result of 
Eq. I|D2|I if the waiting time, r satisfies At — it. In 
Ref. 7J it has been pointed out that a static coupling 
between the resonator and the CPB leads to the desired 
shift of the energy splitting of the charge qubit, 



E,, 



E. 



g,ri 



Ej 






l)Ej 



[huf 



(D4) 



The number state dependent part of the energy differ- 
ence, nA = n4|ApiJj/(£'j — {hvY), can therefore be 
used to implement U . For the parameter values given in 
the previous parts of this paper (A « 10 MHz, v « 100 
MHz, Ej « 10 GHz) the shift A = 40 kHz is actually too 
small but by optimizing the system parameters values for 
A up to 4 MHz pI are possible. With such a system the 
evolution, U can be performed within the decoherence 
time, T2 of the charge qubit |23| which would allow the 
implementation of this fast tomography method. 
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